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ABSTRACT 

In order to study the thermodynamic properties of brane-antibrane systems, 
we compute the finite temperature effective potential of tachyon T in this system 
on the basis of boundary string field theory. At low temperature, the minimum 
of the potential shifts towards T = as the temperature increases. In the D9-D9 
case, the sign of the coefficient of |T| 2 term of the potential changes slightly below 
the Hagedorn temperature. This means that a phase transition occurs near the 
Hagedorn temperature. On the other hand, the coefficient is kept negative in the 
Dp-Dp case with p < 8, and thus a phase transition does not occur. This leads us 
to the conclusion that only a D9-D9 pair and no other (lower dimensional) brane- 
antibrane pairs are created near the Hagedorn temperature. We also discuss a 
phase transition in NS9B-NS9B case as a model of the Hagedorn transition of 
closed strings. 



1 E-mail address: khotta@gauge.scphys.kyoto-u.ac.jp 



1 



1 Introduction 



Recently, significant effort has been devoted to studying non BPS configurations of branes 
such as brane-antibrane pairs and non BPS D-branes (for a review see e.g. [1]). If we 
consider a coincident brane-antibrane pair, the spectrum of open strings on it contains a 
complex tachyon field T. This indicates that the brane-antibrane system is unstable, and 
that the tachyon rolls down from the perturbative vacuum of open strings. If we assume 
that the tachyon potential has a non-trivial minimum, it is expected that the tachyon falls 
into it. Sen conjectured that the negative energy density from the tachyon potential exactly 
cancels the sum of the tension of the brane and the antibrane [2]. This means that the 
brane-antibrane pair disappears and the energy density of this configuration vanishes at the 
potential minimum. In fact, the tachyon potential on this configuration calculated by using 
cubic string field theory approves Sen's conjecture [3]. 

The tachyon potential on the Dp-Dp system in type II string theory was also com- 
puted on the basis of boundary string field theory (BSFT) [4] [5] by Minahan and Zwiebach 
[6] and by Kutasov, Marino and Moore [7], and its qualitative features agree with Sen's 
conjecture. If we denote the complex scalar tachyon field by T, then the potential is given 
by 

\/(T) = 2r p Vexp(-8|T| 2 ), (1.1) 

where V is the volume of the system that we are considering, and t p denotes the tension of 
a Dp-brane, which is defined by 

Tv= To ' (L2) 

where g s is the coupling constant of strings. This potential has the minimum at \T\ = oo. 
We adopt the natural unit c — h — 1 but explicitly keep the slope parameter a'. 

It is interesting to investigate the finite temperature system of Dp-Dp pair. This is 
because there is a possibility that T = becomes a stable minimum at high temperature, in 
analogy with high temperature symmetry restoration in field theories. For this purpose, we 
must compute the finite temperature effective potential of the brane-antibrane system based 
on BSFT, when the complex tachyon field develops its vacuum expectation value. 

Danielsson, Giiijosa and Kruczenski evaluated the finite temperature effective po- 
tential when there is only the tachyon field [8]. The effective action of tachyon field at zero 
temperature can be deduced from the tree level BSFT of open superstrings as [6] [7] 



-1Qt p J d p+1 x 



2a'e- 8 \ T \ 2 dTdT* + -e~^ 



1.3) 



In order to obtain the standard form of the kinetic term, we perform a field redefinition 

1 [ 2T e -\u\* du 
V7T JO 

Then the action (1.3) becomes 



S = Sna'Tp J d p+1 x 
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The mass square of the field is given by 

m2 = w(i e " 8|TW|I )4( |rw|2 -D- (L5) 

and from this the authors of [8] evaluated the finite temperature effective potential in the 
ideal gas limit. It is easy to show from this potential that at sufficiently low temperature the 
potential minimum shifts toward the origin of the complex T plane, T = 0, as the temperature 
increases. 2 On the other hand, computing the potential in the high temperature limit (using 
the standard method [9]), we can show that the critical temperature, at which the potential 
has its minimum at T = 0, is enormously higher than the Hagedorn temperature in the weak 
coupling limit. This is because the tension of the Dp-brane is very large if string coupling is 
very small as we can see from (1.2), and tachyon field must have extremely high energy in 
order to create the Dp-Dp pair. 

There is another approach to the finite temperature Dp-Dp system. Huang at- 
tempted to evaluate the finite temperature effective potential based on BSFT of bosonic 
strings [10]. He concluded that the high temperature tachyon potential has the same form 
with that of the zero temperature up to a temperature dependent string tension, and the 
tachyon rolls down towards the vacuum at T = oo. 3 However, even if we set T = in the 
potential, we cannot reproduce the free energy of on-shell open strings. Therefore, we need 
to compute the finite temperature effective potential in the presence of a constant tachyon 
background based on BSFT and statistical mechanics. Although we will focus on a coincident 
Dp- Dp pair in type II string theory, generalization to non BPS branes is straightforward. 

In addition to the analogy of symmetry restoration of field theories, we have another 
motivation to study the finite temperature Dp-Dp system. The free energy of strings diverges 
above the Hagedorn temperature in general [12]. However, in ten dimensional spacetime, 
we can reach the Hagedorn temperature for closed strings by supplying a finite amount of 
energy, while we need infinite energy to reach the Hagedorn temperature for open strings 
[13]. From this it has been said that, in closed string case, the Hagedorn temperature is 
associated with a phase transition in analogy with the deconfining transition in QCD, while 
it is a 'limiting temperature' in open string case. For this phase transition of closed strings, 
Sathiapalan [14], Kogan [15] and Atick and Witten [16] have argued that the 'winding modes' 
of Euclidean time direction becomes tachyonic above the Hagedorn temperature and a phase 
transition takes place due to the condensation of these tachyon fields. This phase transition 
is called the Hagedorn transition in closed string theory. However, so far we have not known 
where the minimum of this tachyon potential is and in what kind of backgrounds the system 
becomes stable. So we want to shed some light on this problem by investigating the finite 
temperature brane-antibrane system. 

This paper is organized as follows. In §2 we give the free energy of open strings on 
the Dp-Dp system based on BSFT. In §3 we study the finite temperature effective potential 
at low temperature by using the canonical ensemble method. The system near the Hagedorn 
temperature is investigated in §4 and §5 by using the microcanonical ensemble method, 

2 We can show this in a similar way as we will perform in §3. 

3 In a paper [11], he investigated a finite temperature system of parallel brane and antibrane with a finite 
distance and concluded that it becomes stable at sufficiently high temperature. 
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because we cannot use the canonical ensemble method near the Hagedorn temperature as 
we will see in §4. In §6, we discuss a relationship between the Hagedorn transition of closed 
strings and the phase transition on the NS9B-NS9B system. §7 presents our conclusions 
and directions for future work. In appendix A, we describe the finite temperature effective 
potential evaluated by using the canonical ensemble method. 



2 Free Energy of Open Strings 

It is convenient to investigate the finite temperature effective potential by using the canonical 
ensemble method. In this method, the potential of the Dp-Dp system is given by the sum 
of the zero temperature tachyon potential and the free energy of open strings: 



where F is the free energy of open strings. Thus, we begin by considering the free energy of 
open strings in this section. 

We can compute the free energy by using Matsubara formalism. The free energy 
is given by the path integral of connected graphs of strings on the space where Euclidean 
time direction is compactified with the circumference of inverse temperature (5. Let us 
take the weak coupling approximation and treat strings as an ideal gas, that is, we ignore 
the interactions of open strings. We take into account only one-loop amplitude, where the 
corresponding cylinder world-sheets wind the Euclidean time direction at least once. 

There have been attempts to generalize BSFT to the one-loop amplitude of open 
strings on the Dp-Dp system [17] . 4 However, there is an ambiguity in choosing the Weyl 
factors of the two boundaries of a one-loop world-sheet, because the conformal invariance is 
broken by the boundary terms in BSFT. Andreev and Oft have proposed the form of the 
one-loop amplitude of open strings on the Dp-Dp system [18]. They deduced it from the 
principle that its low energy part should coincide with that of the tachyon field which obeys 
the action (1.3). This amplitude can be obtained straightforwardly by choosing boundary 
terms as those of a cylinder world-sheet 



Here we restrict ourselves to the constant tachyon field, which we denote by T. The 1-loop 
amplitude in such background is then given by 



4 Strictly speaking, we must consider quantum master equation in order to investigate quantum corrections 
to the potential. However, the quantum master equation is still missing for BSFT. Thus, we will not discuss 
this subject. 



V(T,P) =2r p Vexp(-8|T| 2 ) + F, 



(2.1) 




(2.2) 




(2.3) 
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Figure 1: For clarity we have displayed the brane and antibrane to be separated. But we 
are considering the coincident brane-antibrane pair. 



We can obtain the same amplitude when we consider the field theory which has the 
mass spectra 

M NS 2 = ^ / (n b + N ns + 2\T\ 2 -'^), (2.4) 

M R 2 = -(N B + N R + 2\T\ 2 ), (2.5) 
a' v 1 

where M^s and Mr are the mass of the Neveu-Schwarz and Ramond sectors, respectively, 
and Nb-i N^s and Nr are the oscillation modes of the boson, Neveu-Schwarz fermion and 
Ramond fermion, respectively. The lowest mode of the NS sector (2.4) coincides with that 
of the tachyon field (1.5). It is convenient to express the free energy by the proper time form 
[19] to compute the free energy from the above mass spectrum. It is given by 

The difference in the signs of each terms come from the fact that the NS sector represents 
spacetime boson and the R sector spacetime fermion. As depicted in the Figure 1, there are 
four types of open strings, namely, two types of the strings whose two ends attach to the 
same brane and two types of the strings which are stretched between the Dp-brane and the 
Dp-brane. If we denote the open strings going from a Dp-brane to a Dp-brane as (Dp, Dp), 
others are denoted as (Dp, Dp), (Dp, Dp) and (Dp, Dp). We must impose the GSO projection 
for (Dp, Dp) and (Dp, Dp), and the opposite GSO projection for (Dp, Dp) and (Dp, Dp). So, 



instead of imposing GSO projection, we may multiply two as an overall factor. Summations 
over the mass can be rewritten by using d function as 



]T exp(-2WM TO 2 r) = 327T 4 e- 4 " |T|2 



M NS A 



^exp(-2 TO 'M R 2 r) = 327r 4 e -^ T l 



^i'(0|ir) 

2 T ($ 2 {U\iT) 



^i'(0|ir) 



Summations over r can be rewritten by using i? function as 



r=l 

oo 

E(-l) r exp 



r 2 /3 2 



r=l 



8na'T / 
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8tkx't 
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Thus, we obtain 5 

F(T,P) : 



16tt 4 V f 00 rfr , , 

tf 3 (0|ir 



(27ra') 



^i'(OKr) 
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' ^(Ojir) 
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?7r 2 a'r 



1 



(2.7) 



(2.J 



(2.9) 



We can reproduce the free energy of on-shell open strings on the Dp-Dp system by substi- 
tuting T = 0. We will discuss the finite temperature effective potential at low temperature 
based on this free energy in the next section. 



3 Finite Temperature Effective Potential at Low Tem- 
perature 

As mentioned in §1, the finite temperature effective potential of the tachyon field model is 
investigated by Danielsson, Giiijosa and Kruczenski [8]. From their calculation we can find 
out that the potential minimum of the finite temperature tachyon potential shifts toward 
T = as the temperature increases. In this section, we will show that we can obtain the 
same result by calculating the finite temperature effective potential of the Dp-Dp system at 
low temperature. 

In the large (3 limit, the integral in the free energy (2.9) can be approximated by the 
large r part. This is because the d function depending on (5 can be expanded in 



8na'T J 

5 We can represent this free energy as the propagator of closed strings by using the modular transformation 
of functions [20]. 
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and these terms are suppressed exponentially unless r is very large. Expanding d functions 
and extracting the leading term in large r region, we obtain 



F(T,ff)~-^- jJVr^exp ^-7r(4|T| 2 - 1 



r 



87ia'r 



where f3 H is the inverse of the Hagedorn temperature 

P„ = 27Tv / 2a 7 . (3.1) 

We should notice that this term corresponds to the contribution of the lowest energy (tachy- 
onic) mode. Above equation can be rewritten by using the third kind of modified Bessel 
function 

£+1 
2 



This free energy coincides with that of the tachyon field model when p — 9 [8] . If we assume 
that both |T| and j3 are very large, the modified Bessel function can be approximated as 

K v {z) ~ Jf z e~\ (3.3) 

and the free energy becomes 

F(T,P) ' exp --f |T| . (3.4) 

We can obtain the finite temperature effective potential by substituting this free energy to 
(2.1). At the potential minimum, this potential satisfies 

and we can obtain the condition for the temperature by substituting the potential we com- 
puted. Taking into consideration the assumption that both \T\ and (5 are very large, this 
gives 

|T| #, (3.6) 

where T denotes the temperature T = and T H denotes the Hagedorn temperature 

T H = ji^ 1 = 1 -=. (3.7) 

2^72^ 



From this we can see that \T\ decreases with increasing T . This implies that the 
potential minimum moves towards \T\ = as the temperature increases, like the tachyon 
field model. 
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4 Complex Temperature Formalism 



It is well-known that perturbative strings have a maximum temperature, which is so called 
the Hagedorn temperature. This comes from the fact that the degeneracy of oscillation 
modes of a single string increases exponentially with the energy, and the density of states 
Q(E) of strings behave as 

fi(£)~e feS , (4.1) 

for large E. The energy of strings is dominated by very high oscillation modes of a single 
high energy string near the Hagedorn temperature. Consequently, we cannot raise the tem- 
perature by adding more energy to the system because a high energy string soaks up most 
of the energy. 

Unfortunately, we cannot use the canonical ensemble method near the Hagedorn 
temperature [21]. We can see this from the relation between density of states and partition 
function, that is, the Laplace transformation 

Z((3)= / dEn(E)e-P E , (4.2) 
Jo 

The integrand is not sharply peeked function of E near the Hagedorn temperature, and the 
canonical ensemble method is not guaranteed to be valid. Thus, we must derive the ther- 
modynamical quantities near the Hagedorn temperature from the microcanonical ensemble 
method, which is more fundamental than the canonical ensemble method in the sense that 
it is derived directly from ergodic theory. This process was performed in the case of closed 
strings [23] and in the case of open strings on D-branes [24]. We will apply this method in 
the case of open strings on the Dp-Dp system. 

We can derive all the statistical variables from the density of states fl(E) in the 
microcanonical ensemble method. Q(E) is given by the inverse Laplace transformation of 

Q(E) = f L+l °° ^Z((3y E , (4.3) 

which is the opposite relation of (4.2). We must choose a constant L such that the path 
of the integral in the complex /5-plane lies on the right side of all the singularities of the 
integrand. Then, we can deform the contour to the left side such as to pick picking up the 
singularities, as is sketched in Figure 2. We will first investigate the singular property of 
the partition function on the complex (3 plane in this section, and then calculate the finite 
temperature effective potential in the next section. 

A general property of the partition function of strings is that it has the leading 
singularity at (5 — f3 H , which we call the Hagedorn singularity. In order to see this, let 
us see the behavior of the free energy near the Hagedorn temperature. In contrast to the 
low temperature case, we must expand the integral in small r region. This can be easily 
performed by making the variable transformation 

1 
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Figure 2: Complex /3 plane. 



and considering large t region. Using the modular transformation of $ functions, we obtain 



F(T,/3) 



16vr 4 V 



p+i 
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y at t 2 exp I - 



fi-Y 



^i'(o|it); 



47r|T| 







^ 4 (0Kt) \ 

AW*)/ 



8vr V / 



87r 2 a' 



(4.4) 



Expanding i? functions and extracting the leading term in large t region near the Hagedorn 
singularity, we obtain 6 



F(T,0) 



4V 



Ph P+l Jo 



p—9 



dt t 2 exp 



47r|T| 



+ 7T 



P 2 



lira' 



t 



(4.5) 



Since we will compute the |T| 2 term of the finite temperature effective potential in the 
vicinity of T = 0, let us expand the free energy in |T| 2 and keep the lower order terms 



F(T,/3) 



4V 



P« p+1 Jo 



E=2 / Ph 2 -P' 



dt t 2 exp 7r 



16ttV|T| 2 /-°° ( 0„ 2 -8 



p+i 



(it t 2 exp 7r- 



A. 



(4.6) 



'Strictly speaking, we must introduce a low energy cutoff as we will see in appendix A. 
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From this we will discuss the finite temperature effective action by using the canonical 
ensemble method in appendix A. If p ^ 9, the integration over t can be easily done by using 
the formula: 

[ f OO 

and we obtain 

r,,m m 4 V f°° , U ' L ^ E 167rV|T| 2 r°° , u^r 

In order to investigate the behavior of the partition function, let us introduce 

W(T,P) = logZ=-pF, (4.9) 

and investigate the singular part of W, which we denote W S i ng . W has a branch point at 
(3 = (3 H , and we can draw a cut to the negative direction of the real axis. Let us evaluate the 
discontinuity of W across the cut. From the definition of the discontinuity, it can be written 

as 

AW = ton [W(T, P + ie)- W(T, [3 - is)} ((3 < (3 H ) 
4V Z" 00 

du u a lim 



P H 7 ~ 2a r(a + 1) Jo e'^+o y u - 2tt« + is' u- 2tt^ 

lQnV\T\ 2 /■«> , ^ / 1 1 



poo 

■ / du u a+l lim 

JO 



(3 H 7 - 2a T(a + 2) Jo e'^+o y u - 2ir£f^ +ie' u- 2rr^f^ - is' 

where e' = 2ne/ ' (3 H and we have introduced 



a = (4.10) 



By using the definition of 5 function 



Six) ee — - — lim ( — — ) (4.11) 



we can rewrite AW as 



r(a + l)/i/ H r(a + 2)^ 

This non-zero AW guarantees that /3 = /3 H is a branch point. We can deduce W sing from 
AW as follows [23] [24]. 



(a) a : half-integer (p : even) 



(P — (3 H ) a has discontinuity at P — ft H when a is a half-integer. Using this fact and a 
formula for V function 

r(z)r(i - z) = (4.13) 

sin irz 

10 



we can express W S i ng as 

W sing ~ [C a (P - (3 H ) a - D a \T\ 2 ((3 - (3 H ) a+l ]V, (4.14) 

where 

C = (4,5, 

D . = ^-') . (4,0) 
(b) a : integer (a ^ — 1) (p : odd (p ^ 9)) 

In this case, (f3 — f3 H ) a has no discontinuity. However, \og[f3 H ~ l (f3 — (3 H )] has disconti- 
nuity at (3 — f3 H and we can express W S i ng as 

W smg ~ [C Q (/3 - 0*)° - L» Q |T| 2 (/3 - /5 H ) a+1 ]V log f^^ 1 ) , (4-17) 



where 



(-l) a+1 4(2 7r) a 

r(a + 1)& 



^ = / t,,- (4.18) 



£) = (- 1 ) Q+28 ^) Q+2 (419) 



(c) a = -1 (p = 9) 



In this case, we can perform the integral in the first term of (4.6) explicitly. For the 

second term, we can obtain the discontinuity in the same way when a is an integer. 
Combining two terms, we can express W sing as 

^~^- D -W>cg(^), (4.20) 

where 

°- = ^7- (4 ' 21) 

= -£5- (4.22) 

Ph 

The first term gives the pole in the complex /3-plane. 



We thus have obtained W sing for all possible a, which can summarized as 

W sing ~ C a V(/3 - /k)°(l°g (^p)) -^V|T| 2 (/3-^r (^log^^^^ , (4.23) 

where a is 1 if a is an integer except —1 and otherwise, and 6 is 1 if a is an integer and 
otherwise. The singular part of the partition function can be obtained from (4.9) and this 
W sing . We can compute the density of states from it. We will see it in the next section. 
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5 Finite Temperature Effective Potential near the Hage- 
dorn Temperature 



We can derive thermodynamical variables from the density of states in the microcanonical 
ensemble method, as we mentioned in the previous section. The density of states Q(E) can 
be obtained from the partition function Z(/3) by using the inverse Laplace transformation 
(4.3). Z(P) is related to W{(3) as 



Z(T,0) =e w = ex P [W reg + W 1 



smg\ ■ 



(5.1) 



Here W reg is the regular part of W, which can be expanded in a power of (f3 — f3 H ) as 



reg 



\ V-a V(P-P H ) + O{V(P-P L 



(5.2) 



It is easily to see that Ao is the constant which has the dimension of the inverse of volume, 
and do is the constant which has the dimension of an energy density. Thus, using eqs. (4.3), 
(4.23), (5.1) and (5.2), we obtain the density of states as 



Q(T, E) 



dp 
c a 27ri 



exp 



(P-f3 H )E + 0(V([3-[3 H ) 2 ) 

'P-p H 



+C a V(P-p H ) a log 



Ph 



-D a V\T\ 2 (P-P H ) a log 



(5.3) 



where E = E — a V. Furthermore, knowing the explicit form of Q(E), we can calculate the 
entropy S(E) as 

S(T, E) = log Q(T, E)8E, (5.4) 
and the inverse temperature P as 

"=!• _ ^ 

Then, the finite temperature effective potential of the Dp-Dp system is given by 



V(T,E) = 2r p Vexp(-8|T| 2 ) - p~ x S. 



(5.6) 



Now, we are ready to calculate the finite temperature effective potential by using the micro- 
canonical ensemble method. We will give a discussion for each value of p separately, since 
our approximation to calculate Cl(T, E) differs for each value of p. 

(a) p = 9 (a = -1) 



Let us first consider the case of the D9-D9 system, which is the most interesting case. 
In this case (5.3) is rewritten as 



n(T,E) ~ e 



dP_ ( P-Pn 

c a 2m I p H 



-D_iV|T| 2 



exp 



(P-P H )E + 



c_y_ 



(5.7) 
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Let us suppose that both E and EjV are very large. Then, the saddle point method 
works well because the exponent in the integrand is very large. The result is 7 



n(T,E) ~ 



2v^F \c_y 



exp /3 H £ + A V + 2VC_ 1 V£7 . (5.8) 



In (5.7), we ignored 0(V((3 — I3 H ) 2 ) term when we derive this equation from (5.3). We 
will examine the correction of the saddle point method from this term later in the 
p = 6 case. Substituting Q(E) into (5.4), we obtain 



-lo ( E 



+(3 H E + A V + 2JC^VE, 



and from (5.5), we get 



D_y\T\ 2 3 /C^V 
4l + ^ + V^- 



(5.9) 



(5.10) 



The finite temperature effective potential can be derived from (5.6). In order to argue 
the stability of the brane-antibrane system, we need only |T| 2 term of V(T,E). This 
term is given by 

8ttV , fn(3 H w W 



-16r 9 V + 



log 



T 



(5.11) 



f3 H w '~°\ 2V 

It should be noted that the second term in the coefficient of |T| 2 increases with in- 
creasing E. Since the first term is constant as far as V and r 9 fixed, the sign of the 
\T\ 2 term changes from negative to positive at large E. The coefficient vanishes when 



E 2V 



'2(3 H 10 r 9 



7T(3 H 



10 



cxp 



7T 



If we approximate (5.10) as 



(5.12) 



(5.13) 



we can derive the critical temperature T c at which the coefficient vanishes. The result 
is 

Ph Tg 



1 + exp 



7T 



(5.14) 



Here we see that this temperature is very close to the Hagedorn temperature since r 9 is 
very large if the coupling of strings is very small. Above this temperature, the coefficient 
of |T| 2 is positive and T = becomes the potential minimum. This implies that a 
phase transition occurs at the temperature T c which is slightly below the Hagedorn 
temperature, and the D9-D9 system is stable above this temperature. 



7 We can obtain the same result by using modified Bessel function. 
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p= 8 ( a = -1/2) 



Next, let us consider the case of the D8-D8 system, where we will see that a result 
strikingly different from the case of the D9-D9 system arises. From (5.3) the density 
of states is given by 



Q(T,E) ~ e 



I3 H E+\ V 



/ — exp 

Jc 



o P H E+\ V 



iCa 2ni 
dp_ 

c a 2-Ki 



(P - (3 H )E + C_i V{(3 - (3 H )~ - D_,V\T\\P - 

l-D_ h V\T\ 2 (P-p H y\ 

(J3 - i3 H )E + C_ h V(J3 - j3 H )-t] , (5.15) 



2 

x exp 



where we take the small |T| 2 approximation in the second equality. We can also use 
the saddle point method and obtain 



n(r, e) 



c_iM 
2 



exp 



P H E + A V + 



ZC.iWe* C_i^D_iVhT\ 2 



V 



2§ 



. (5.16) 



We can calculate the entropy S, the inverse temperature f3 and the potential V(T, E) 
from (5.4), (5.5) and (5.6) as in the D9-D9 case. The |T| 2 term of V(T, E) is given by 



-16r 8 V - 



2 fvf 



(5.17) 



3p H LZ E' 

It should be noted that the second term in the coefficient of |T| 2 decreases as E gets 
large. Thus, the coefficient of |T| 2 remains negative for large E. This implies that a 
phase transition does not occur unlike in the D9-D9 case. 



p = 7 (a = 0) 



In this case, the density of states can be obtained from (5.3) as 

D /3 H E+\ v f d@ 



n(r, e) ~ e 



c a 27ri 



: exp 



(P - p H )E + (C - A)m 2 (/3 - P H )) Vlog 



' P ~ P, 



(5.18) 

and we can also use the saddle point method if we expand in |T| 2 like in the D8-D8 
case. The result is 



Q(T,£)^^^exp 



lixE 



P H E + A V - C V + (C V + CoD °l |T| j log (-g 



(5. 



9) 



We can calculate the finite temperature effective potential from this, and its |T| 2 term 
is given by 

2 7 tt 2 V 2 , fp H 8 W 



-16r 7 V 



Ph 16 E 



4V 



\T\ 



(5.20) 



From this we can see that the coefficient is always negative, so that a phase transition 
does not occur. 
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p= 6 (a = 1/2) 



In this case, the density of states can be obtained from (5.3) as 
j h e+x v f d P 



Q(T, E) ~ & 



(5.21) 

and we can also use the saddle point method if we expand in \T\ 2 like in the D8-D8 
case. The result is 



tt(T,E) 



CiV 

2 



exp 



(3 H E + A V 



Ci 2 V 2 CilDiV 4 
~w + 



4£ 



8E 



(5.22) 



We can calculate the finite temperature effective potential from this, and its |T| 2 term 
is given by 



■16r 6 V + 



2 17 tt 6 V 4 
{3 H 28 E 3 j 



ITI 



(5.23) 



From this we can see that the coefficient is always negative, so that a phase transition 
does not occur. 

Here, we must consider corrections to the saddle point approximation. Such corrections 
come from 0(V(P — I3 H ) 2 ) term in (5.3), which we have ignored so far. This term gives 
rise to a shift of the saddle point. Let us consider in the case that a is a half integer. 
If we write this term as 



n=0 

then the correction enters as a multiplicative factor 



cxp 



£0 V 



n+2 \ 
a-1 



n=0 



V 



(5.24) 



(5.25) 



If the ratio E/V is very large, we can use the saddle point method for a < 1/2, 
because the power of E/V is negative and the correction becomes negligible. However, 
we cannot ignore the correction for a > 3/2, because the power is positive and the 
correction becomes very large for sufficiently large n. Similar calculation leads us to 
the result that we can use the saddle point method also in the case of a — 0, — 1. If 
a is a positive integer then 0(V(f3 — (3 H ) 2 ) term is of the same order to other terms. 
Therefore, we cannot use the saddle point method when a > 1. 

p = 5 (a = 1) 

In this case, we cannot use the saddle point method, and thus we adapt another 
approximation as follows. The density of states can be obtained from (5.3) as 



n(T,E) 



o l3 H E+X V 



r dp 
Jc a 2m 



(P-Ph)e 



(d(/3 - p H ) - D X \T\ 2 {(3 -(3 H ) 2 )V log 



(5.26) 
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If we define 

then (5.26) can be rewritten as 

e P H E+\ n V 



(5.27) 



n(r,E) ~ 



E' 



r dz 
Jc b 2ni 

x exp 



dVz 



E> 



log 



D x V\T \ l z 



Vz 

2„2 / 



log 



Til 



Til 



where we define 



E' = E-C X V log 



V 



5 ' 



(5.28) 



(5.29) 



for simplicity. The contour of the integral is deformed in the complex plane, as sketched 
in Figure 3. We approximate the integral in (5.28) by the contour integral in both side 
of the cut as 



r dz 
Jzi 2ni 

I 

r z i dz 

JO 7T 



exp 



C x Vz 



f Zl dz -z 

' 1 2~Ki' 6XP 



E' \ 
C{Vz 



log 



exp 



x sin 



E' 

CxVz. 

ioe 

E' & 



log 



_Vz_ 

Ph'e*, 

Vz 



— m — 



DMTVz 



2^2 



3hE< 



m — 



E' 2 
D l V\T\ 2 z 2 



log 



— m 



Vz 



DiV|T| V 
E' 2 



p H 6 E'J 

ixCxVz nD 1 V\T\ 2 z 2 ^ 



E' 2 
log 



log 



+ 7ri 



E> 



+ 



E' 2 



where 

Zl = -(L'-p H )E'. (5.30) 

Since we are assuming that E is very large, the integrand can be expanded in l/E', 
and we can take a limit zi ^ oo. The integral becomes 



C X 2 V 2 , 

9— l°g 

E' 



V 



DiVm 



A r E', 



E> 



r(3) + ^Pr'(s), (5.3i) 



where we express the derivative of V function as T'(z). If we include 0(V([3 — f3 H ) 2 ) 
term in (5.3), it can be rewritten as 0(Vz 2 / E' 2 ) and is a function of ((3 — /3 H ) to the 
integer power, so that we can drop it in this expansion. We can calculate T'(z) from 
the relation between T function and digamma function ip(z) 



1>(z) 



r\z) 

T(z) 



and from 



m = 2-7, 



(5.32) 



(5.33) 



16 



c s Li. 



* r > 



o 



Figure 3: Complex z plane. 



where 7 is the Euler constant 7 = 0.57721 • --. Substituting the result into (5.28), we 
obtain 



MT,E)~-?^e?" E+ *> v 
V ' p H 6 E' 2 



1 - 



16vrV 



7 log 





8tt 2 


T 


{p H 6 E'J 


(3 H E> 



Ph E' 



^ (5.34) 

We can calculate the finite temperature effective potential from this, and its \T\ 2 term 
is given by 

16vr 2 

-16t 5 V + - 77 ^ 



\T\ 



(5.35) 

From this we can see that the coefficient is always negative, so that a phase transition 
does not occur. 

p = 4,2,0 (a = 3/2,5/2,7/2) 



We can calculate the finite temperature effective potential in analogy with the p = 5 
case. The density of states can be obtained from (5.3) as 



n(T, E) ~ e 



/3hE+X V 



dp 



ICa 2ni 
Using the transformation 



exp 



(P - P„)E + (c a (p - p H y - D a \T\\p - p H ) a+1 ) V 

(5.36) 



<P-p H )E, 



(5.37) 
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we get 

n(r,E) 



e (3 H E+\oV , dz 



E 



r dz 
Jc b 2iri 



exp 



H -7 a D lTl 2 r Q + 1 \ 



-=rv 6 



E 



E 



■a+l 



(5.38) 

The contour of the integral is deformed as in Figure 3, but Z\ is here replaced by z 2 , 
where 

z 2 = -(L'-p H )E. (5.39) 
We can approximate the integral by the contour integral in both sides of the cut as 



« dz _ z . (C a sm(-7ra)Vz a _ D a sm(-7i(a + l))V\T\ 2 z 
e sin I =jj 

7T \ E 



2„a+l ' 



E a+l 



(5.40) 



Expanding the integrand in 1/E and taking the limit z 2 — > oo, the integral becomes 
4(2vr) Q r(-a)r(a + 1) sin(-Tra) V 



7-a 



j(27r)"+ 2 r(-a - l)r(a + 2) sin(-7r(q + 1))|T| 2 V 



Using (4.13), we obtain 

fi(T, £) 



(5.41) 



(5.42) 



We can calculate the finite temperature effective potential from this, and its |T| 2 term 
is given by 

16tt 2 



I T I ^ 



(5.43) 



From this we can see that the coefficient is always negative, so that a phase transition 
does not occur. 



p = 3,l (a = 2,3) 



We can also calculate the finite temperature effective potential in a similar way as that 
we have done in the above two cases. The density of states can be obtained from (5.3) 

as 



tt(T,E) 



o /3 H E+X V 



dp 
c a 2ni 



- cxp 



(P - Ph)e 



+ (C Q (/3 - p H ) a - D a \T\ 2 (P - p H ) a+1 ) V log 
and a similar calculation leads us to the same result, that is, 



' P ~ & 



(5.44) 



Ph E \ 



8n 2 \T\ 2 ' 
PhE 



(5.45) 
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We can calculate the finite temperature effective potential from this, and its |T| 2 term 
is given by 

T| 2 . (5.46) 

From this we can see that the coefficient is always negative, so that a phase transition 
does not occur. 

Let us summarize the results of this section. In D9-D9 case, the sign of coefficient of 
the |T| 2 term of the finite temperature effective potential changes from negative to positive 
near the Hagedorn temperature as the temperature increases, while it remains negative in the 
case of Dp-Dp with p < 8. These results lead us to the conclusion that a phase transition 
takes place only in the case of D9-D9 pair. Therefore, not the lower dimensional brane- 
antibrane pair but the space-filling D9-D9 pair is created near the Hagedorn temperature. 



6 Speculation about the Hagedorn Transition 

Let us turn to discuss the Hagedorn transition in closed string theory. If we consider the 
closed strings near the Hagedorn temperature, the strings are highly excited and the energy 
density is very large. Therefore, the effective coupling of strings becomes very large. As a 
consequence, we must consider this transition in the strong coupling region. 

In addition, it seems correct that the coupling of strings becomes very large during 
this phase transition for the following reason. Let us recall the relation between type II 
theory and type theory [25]. The one- loop free energy of type II theory in the high tem- 
perature limit coincides with the one-loop amplitude function of type theory. On the other 
hand, Bergman and Gaberdiel conjectured that type OA theory is obtained by the orbifold 
compactification of M-theory [26]. The coupling constant of type OA theory becomes large as 
the radius of the compactified eleventh dimension increases. While type OA theory includes 
a tachyon field and is unstable, M-theory has no tachyon field and is stable. The above 
conjecture implies that type OA tachyon becomes massive at sufficiently strong coupling, 
and that the tachyonic instability disappears in the strong coupling regime. Therefore, the 
weak coupling type OA theory undergoes a phase transition by tachyon condensation and the 
eleventh dimension decompactifies completely. From this viewpoint one may say that the 
coupling constant increases unlimitedly during the Hagedorn transition in type IIA theory. 
For type IIB theory, similar speculation can be deduced from the relation between type OB 
theory and M-theory. Thus, we must discuss the phase transition in the strong coupling. 

Let us consider the phase transition in type IIB theory. This theory has the S-duality 
symmetry which interchanges the strong coupling region with the weak coupling one. In the 
previous section, we have shown that the D9-D9 system becomes stable near the Hagedorn 
temperature. If we apply the S-duality to D9-brane, we obtain an NS9B-brane [27]. The 
world volume theory of D9-brane is described by open fundamental strings, which live on 
it. Since the S-dual object of fundamental string is D-string, the world volume theory of 



16tt 2 
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NS9B-brane can be described by the open D-strings. D-strings can be treated perturbatively 
in large g s region, where g s is the coupling constant for closed fundamental strings. 



Let us see the Hagedorn temperature for D-strings. The Hagedorn temperature for 
fundamental strings depends only on their tension, which is given by 

T - = db- (6 - 1) 

On the other hand, the tension of the D-string can be obtained by substituting p = 1 into 
(1.2) as 8 

Tm = 7T^- ( 6 - 2 ) 

Therefore, we can obtain the Hagedorn temperature for D-strings by replacing a' in the 
Hagedorn temperature (3.7) with a'g s . The result is 

T DH = 1 . (6.3) 
2iTy/2a'g s 

We can see that the Hagedorn temperature for D-strings is very small when g s is very large. 
Our calculation implies that the critical temperature for the phase transition in the case of 
NS9B-NS9B is slightly below this temperature and it is very small. The reason for this can 
be understood from the tension of the NS9-brane, which is given by [27] 

r ns9b = — — r— — -f — j. (6.4) 
(2vr)VV 



From this we can see that the NS9B-NS9B pair is created in the low energy region when g s 
is very large. 

From these observations we propose the following model of the Hagedorn transition. 
Let us start from closed fundamental strings at low temperature. Closed strings are excited 
as the temperature increases, and the density of strings become very large near the Hagedorn 
temperature. This lead to the increase of g s , and the tension of NS9B-brane becomes small. 
Then a phase transition occurs, and the NS9B-NS9B system appears. 

Let us make some comments on this model. Firstly, it seems reasonable to assume 
that higher dimensional NS-charged objects are realized as a configuration of infinitely many 
lower dimensional NS-charged objects in analogy with the D-brane case [28]. In this sense, 
one may say that the NS9B-NS9B system consists of fundamental strings in our model. 

Secondly, if we consider the Hagedorn transition in flat spacetime, there is no pre- 
ferred choice of directions. Thus, it is natural to think that closed fundamental strings 
change to the spacetime-filling brane, since the spacetime-filling branes manifestly preserve 
all spacetime symmetries, while the lower- dimensional branes break some of the symmetries 
Our results seem very reasonable in this sense. 



*Since D-string is a BPS object, we may extend this formula to the strong coupling region. 
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Thirdly, the coupling constant of strings is related to the expectation value of dilaton 
field. If we denote the expectation value as <fi then the coupling constant is represented as 

9s = e*. (6.5) 

Since the dilaton field is one of the massless modes of closed strings, the coupling constant is 
determined from closed string side. As a consequence, open D-string are kept weakly coupled 
as long as the closed fundamental string coupling g s is very large. 



7 Conclusion and Discussion 



In this paper, we have discussed the behavior of the finite temperature effective potential on 
the brane-antibrane pair in the constant tachyon background. We evaluated the potential at 
low temperature by using the canonical ensemble method and concluded that the potential 
minimum shifts towards T = as the temperature increases in the low temperature region. 
Near the Hagedorn temperature, we have calculated the \T\ 2 term of the potential by using 
the microcanonical ensemble method. For the D9-D9 system, the sign of the coefficient of 
this term changes from negative to positive at slightly below the Hagedorn temperature. This 
implies that a phase transition occurs at this temperature and the D9-D9 system becomes 
stable above this temperature. This result is in sharp contrast to lower dimensional brane- 
antibrane cases. For the Dp-Dp system with p < 8, the coefficient remains negative near the 
Hagedorn temperature, so that such a phase transition does not occur. We thus concluded 
that not a lower dimensional brane-antibrane pair but a D9-D9 pair is created near the 
Hagedorn temperature. 

This difference comes from the contribution from the momentum modes of open 
strings, which can be taken in p-dimensional spatial directions on Dp-brane. If we compactify 
the background spacetime, we must take into consideration not only momentum modes but 
also winding modes. It would be interesting to investigate the finite temperature effective 
potential in a compact background spacetime because the energy spectrum of momentum 
and winding modes depend on the background spacetime. In particular, we might be able 
to analyze the finite temperature effective action explicitly in the case of torus background, 
which is the simplest non-trivial one. 

Our calculation is based on the one-loop amplitude of open strings, which has been 
proposed by Andreev and Oft [18]. As we have mentioned in §2, we have a problem in 
choosing the Weyl factors in two boundary terms of one- loop world-sheet. Although the 
choice of Andreev and Oft has the meaning that both side of cylinder world-sheet is treated 
on an equal footing, we must explain the fundamental reason why this choice is successful 
even if it is correct. However, we only need to analyze the vicinity of T = in order to 
investigate whether the phase transition occurs or not. Our calculation is valid if the mass 
square is shifted as equations (2.4) and (2.5) for small \T\. 

We have proposed that the Hagedorn transition in type IIB string theory might 
be a transition from closed strings to the NS9B-NS9B system. This proposal comes from 
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our calculation and the S-duality. Since the tension of NS9B-brane decreases as the string 
coupling constant increases, it is natural to think that not only one pair but also large number 
of brane-antibrane pairs appear in this phase transition. Thus we must consider the finite 
temperature system of n NS9B-branes and n NS9B-branes, or their S-dual objects, namely, 
those of n D9-branes and n D9-branes. 

Type IIA string theory is obtained by applying a T-duality to type IIB theory. Since 
the T-dual object of NS9B-brane is NS9A-brane [27], our proposal is replaced by the phase 
transition from closed strings to the NS9A-NS9A system in type IIA theory. In strong 
coupling region, type IIA theory is lifted to eleven dimensional M-theory, and NS9A-brane 
to M9-brane. Thus, we might have to study the properties of M9-brane in order to investigate 
the Hagedorn transition in type IIA theory. 

If there exist charged objects such as D-branes and NS-branes before the Hagedorn 
transition occurs, then it seems natural that NS9-NS9 system includes these objects. Topo- 
logical configuration of these objects is described by (the S-dual of) the K-theory [29] [30]. 
It would be interesting to investigate our model in terms of the K-theory. 

The application of our model to cosmology is very important because brane-antibrane 
pairs might be excited if the temperature of the early universe is close to the Hagedorn tem- 
perature. The tension of the brane-antibrane system behaves as the cosmological constant, 
so that the universe experiences inflation. In fact, much work has been done in this direction 
[31] [32]. 

Finally, the phase transition to the 9-brane is reminiscent of the Plank solid model 
of Schwarzschild black holes [33]. It might be interesting to study the black holes as the 
NS9-NS9 system. 
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A Finite Temperature Effective Potential in Canonical 
Ensemble Method 

As we commented in §4, the calculation based on the canonical ensemble cannot be trusted if 
the temperature is closed to the Hagedorn temperature. In this appendix, we show that the 
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canonical ensemble method gives different results from those obtained by the microcanon- 
ical ensemble method, by explicitly calculating the potential based only on the canonical 
ensemble. From (4.6) the free energy is given by 



4V p-9 ( p H 2 - p 2 

16ttV|T| 2 r ^ / P H 2 ~P 2 



+ a p+i /, dtt— exp(7r ^ f- t), (A.l) 



where we have introduced the low energy cutoff. This cutoff is required when we made an 
approximation to derive (4.5). Using the incomplete T function, (A.l) can be rewritten as 

7T— p H s {p 2 -f3 H 2 )— \ 2 Ph J 

16ttV|T| 2 /p - 9 /3 2 -pV \ 

The potential for each p can be derived from this as follows. 



(a) p = 9 



If we substitute p = 9 into the first term of (A.2), the first argument of incomplete V 
function becomes one, so that we can set A = 0. For the second term of (A.2) we can 
use the following formula for the incomplete T function; 



00 C_ \\ n x n 
r(0,a;) = -7 - log a; - ^ 



n=l 



n ■ ni 



(A.3) 



where 7 is the Euler constant. Combining these two terms, we get 

4V 



F(T,(3) 



" 10 7T — A 



pV 



Ph 



(A.4) 



We can obtain the finite temperature effective potential from (2.1), and its \T\ 2 term 
is given by 



-16r 9 V 



16ttV 



log vr 



P 2 ~ P, 



-A 



I T I ^ 



(A.5) 



P H W S V Pn 

This has a similar form to (5.11) which was derived by using the microcanonical en- 
semble method. The coefficient vanishes at the critical temperature T c given by 



T c = p- 1 ~ 



P Z H + 
Pu- 1 



Ph 2 




7tA 


exp 




1 


1 + 




2ttA 



7T 



exp 



7T 



(A.6) 
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If we set A = (2ir) 1 , the result is the same as (5.14) which is derived by using the 
microcanonical method. 9 



(b) p : even 



In this case we may set A = and we obtain T functions. T function with half integer 
argument is given by 

Substituting this into (A. 2), we obtain 

(-l)5-l2 10 ^(4-l)! (/ 32_ /3ff 2 ) i ?v 



F(T,f3) ~ 



ttI- 4 (8-p)! 



T\ 2 . (A.8) 



We can obtain the finite temperature effective potential from (2.1), and its |T| 2 term 
is given by 



T\\ (A.9) 



(-l)H2^5-f!(^-^)¥v 

10T V V + ■ 



The coefficient remains negative near the Hagedorn temperature, namely /3 ~ f3 H . 
Although the result agrees with (5.17) and (5.25) for p = 8 and p — 6, respectively, up 
to numerical coefficient, we obtain different results for p = 4, 2 , 0. 

(c) p : odd (p ^ 9) 

The incomplete Y function whose first argument is a negative integer can be expanded 

as 

T(-n, x) = -e~ x Ti-iy-'in - s)\ x~ n+s - 1 + ^^r(0, x). (A.10) 
nl ~± nl 

In our case, A — > corresponds to x — > 0, so that incomplete T function can be 
approximated by the last term. Thus, the free energy can be approximated as 

i6^(-i)^ p.* - ^)^v |r|2log L^PlA . (A . n) 



9 In order to keep the argument of the logarithm in (4.20) dimensionless, we have divided {(3 — (3 H ) by (3 H . 
But this denominator is ambiguous as long as it has the dimension of the inverse of the temperature. For 
the different denominator, we can obtain the same result by shifting the cutoff A in the canonical ensemble 
method. 
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We can obtain the finite temperature effective potential from (2.1), and its |T| 2 term 
is given by 



-16r p V + 



10 



log 71- 



-A 



ITI 



(A.12) 



The coefficient remains negative near the Hagedorn temperature either. The result 
agrees with (5.20) for p = 7 if we set A = (2tt)~ 1 . However, we obtain different results 
for p = 5, 3 , 1. 
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